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HE: 异 叶 泽 兰 属于 菊 科 泽 兰 属 ， 是 该 属 分 布 海拔 相对 较 高 的 植物 , 分 布 在 青藏 高 原 东 部 和 
横断 山海 拔 1 700~3 000 m 左右 的 地 区 。 本 研究 利用 ycf6-psbM 和 rp132-trnLl 两 个 叶绿体 DNA 
(cpDNA) 片段 以 及 核 DNA 片段 ITS OITS) 作为 分 子 标记 ,研究 了 异 叶 泽 兰 的 遗传 多 样 
性 及 其 分 布 特征 , 同时 探讨 了 其 居 群 历史 动态 。 叶绿体 片段 联合 分 析 结 果 显示 单 倍 型 多 态 性 
指数 Hd WN 0.656, 核 苷 酸 多 态 性 (rr) 尖 .001 61; 而 TS 的 平 基因 型 多 态 性 指数 Hd 为 0.687, 
PKA IRE APEA 0.002 35。 因 此 异 叶 泽 兰 在 物种 水 平 遗传 多 样 性 水 平 不 高 。cpDNA 和 nITS 
分 析 结 果 都 显示 异 叶 泽 兰 居 群 水 平 总 的 遗传 多 样 性 大 于 居 群 内 平均 遗传 多 样 性 , 遗传 变异 主 
要 发 生 在 居 群 间 , 居 群 间 存 在 明显 的 遗传 分 化 (cpDNA: Gst= 0.679, Nst= 0.655, Fsr= 0.655; 
nITS: Gst = 0.543, Nst = 0.370, Fsr= 0.584)。 但 是 ， 由 于 Nst 值 小 于 Gst 值 ， 异 叶 译 兰 的 
分 布 不 具有 明显 的 谱系 地 理 结构 。 基 于 单 倍 型 地 理 分 布 以 及 Network 分 析 推 测 横 断 山 区 南部 
《川西 南 - 滇 西北 ) 和 云南 中 部 可 能 是 异 叶 泽 兰 在 第 四 纪 冰 期 时 的 两 个 避难 所 ， 中 性 检验 和 
失 配 分 析 的 结果 支持 异 叶 泽 兰 在 冰期 后 未 发 生 过 居 群 扩张 。 
关键 词 ， 异 叶 泽 兰 ， 遗 传 多 样 性 ， 居 群 历史 动态 ， 横 断 山 区 ， 和 避难 所 
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Abstract: Eupatorium heterophyllum, belonging to Eupatorium (Asteraceae), is distributed in the 
eastern Qinghai-Tibet Plateau and Hengduan Mountain regions, and has relatively high elevation 
in comparison with those of other species within this genus. A survey of two chloroplast DNA 
(cpDNA) fragments (ycf6-psbM and rpl32-trnL), and one nuclear DNA fragment ITS (nITS) was 


carried out to assess the genetic diversity and infer the population demography history of E. 
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heterophyllum. The results showed that the haplotype diversity (Hd) of combined cpDNA 
fragments was 0.656, and the nucleotide diversity (z) was 0.001 61 at the species level. These two 
index values of ITS were 0.687 and 0.002 35, respectively. All these data showed that E. 
heterophyllum had relatively low level of genetic diversity. Both cpDNA and ITS data indicated 
that the total genetic diversity of E. heterophyllum was higher than the average value of 
populations, and the genetic variations occurred mainly among populations. Significant genetic 
differentiation exists among populations (cpDNA: Gst = 0.679, Nst = 0.655, Fsr= 0.655; nITS: 
Gst = 0.543, Nst = 0.370, Fsr= 0.584). However, there is no obvious phylogeographical structure 
occurred in E. heterophyllum (Nst < Gst). Both the pattern of haplotype distribution and Network 
structure of haplotypes indicated that the southern Hengduan Mountian regions and central 
Yunnan province probably were the two refugia for E. heterophyllum during the Quaternary 
Glaciation period. And E. heterophyllum has not undergone population expansion after glacial 
period which was supported by the neutral test and mismatch analysis. 
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FEM Y= (Eupatorium heterophyllum) 属于 菊 科 译 兰 属 (Eupatorium), 该 属 曾 被 认为 是 
分 布 在 美洲 及 欧 亚 大 陆 的 菊 科 大 属 , 2J 600 种 , 被 国人 所 熟知 的 外 来 入 侵 物 种 紫荆 泽 兰 就 曾 
被 置 于 该 属 ， 拉 丁 名 为 E. adenophora ( 林 容 等 ，1985; Chen et al., 2011). 但 后 来 该 属 被 重新 
界定 为 仅 分 布 在 欧 亚 大 陆 的 北极 -第 三 纪 分 布 小 属 ， 仅 包含 45 个 种 左右 (King & Robinson, 
1970)， 而 分 布 于 中 南美 洲 的 类 群 不 再 被 归 入 该 属 内 ， 如 紫荆 泽 兰 被 归并 到 紫 茎 泽 兰 属 
Ageratinaq， 定 名 为 Ageratina adenophora (King & Robinson, 1970; Chen et al., 2011)。 被 重 
新 定义 的 泽 兰 属 属 北 温带 类 群 ， 北 美 分 布 有 27 种， 东亚 有 25 种 ，1 种 在 欧洲 ， 是 个 典型 的 
东亚 北美 间断 分 布 类 群 (Schilling et al, 1999; Ito et al., 2000)， 在 中 国 大 陆 和 台湾 地 区 约 分 布 
有 14 种， 其 中 6 种 为 特有 种 (Chen et al., 2011)。 系 统 发 育 和 植物 地 理学 分 析 表 明 ， 该 类 群 
为 北美 起 源 和 分 化 ， 其 一 分 支 大 约 在 6.18~11.6 Ma MEZA) 时 经 白 伶 海 峡 陆 桥 扩散 至 欧 
亚 大 陆 ， 随 后 在 亚洲 经 历 了 辐射 演化 (Schilling et al., 1999; Schmidt & Schilling, 2000)。 其 
中 ， 东 亚 分 布 的 种 类 在 形态 表 型 以 及 染色 体 核 型 上 差异 都 比较 小 (Watanabe et al., 1990), 
物种 间 可 能 存在 自然 杂交 的 现象 (Schmidt & Schilling, 2000)。 在 中 国 分 布 的 约 14 个 种 中 ， 
异 叶 泽 兰 分 布 在 青藏 高 原 东 部 和 横断 山海 拔 1 700~3 000 m 左右 的 地 区 ， 是 泽 兰 属 唯 一 一 个 
分 布 海拔 较 高 的 物种 ， 生 于 山坡 林 下 、 林 缘 、 草 地 及 河谷 〈Chen et al., 2011)， 其 化 学 成 分 
在 不 同居 群 间 存 在 差异 (Saito et al., 2014). 
横断 山区 是 泛 北极 植物 区 系 的 一 个 区 〈 李 锡 文 和 李 捷 ，1993)， 该 地 区 的 种 子 植物 区 系 
基本 上 是 温带 性 质 的 , 尤其 是 北 温带 成 分 ， 这 其 中 许多 又 是 北极 -第 三 纪 植 物 ， 这 类 植物 在 
中 新 世 以 后 伴随 气温 的 急剧 下 降 而 南 迁 。 秦 岭 - 黄 河 一 线 并 环 沿 四 川 倪 地 可 能 是 北极 第 三 纪 
成 分 向 喜马拉雅 -横断 山 迁 移 的 主要 路 线 〈 孙 航 , 2002)。 该 地 区 在 第 四 纪 冰 期 时 未 被 统一 的 
大 冰 盖 覆盖 , 同时 又 受 东 亚 季 风气 候 和 印度 季风 气候 的 影响 , 这 给 许多 植物 类 群 提供 了 避难 
所 ， 尤 其 是 海拔 相对 较 低 的 地 区 (Qiu et al., 2011; 更 吉 卓 玛 等 ，2018)。 在 冰 后 期 ， 横 断 山 
区 成 为 现代 温带 植物 的 重要 起 源 地 和 辐射 地 ( 李 锡 文 和 李 捷 ,1993; 于 海 彬 和 张 镜 锂 ,2013 )， 
物种 多 样 性 和 遗传 多 样 性 都 非常 高 (Yu et al., 2019)。 异 叶 泽 兰 作为 泽 兰 属 在 高 海拔 地 区 分 
布 物种 , 在 第 四 纪 冰 期 及 冰期 后 存在 怎样 的 居 群 历史 动态 , 现 有 分 布 居 群 又 具有 怎样 的 遗传 
多 样 性 分 布 式样 。 本 文 基 于 两 个 叶绿体 DNA 片段 (ycf6-psbM 和 rpl32-trnL) 和 核 DNA Fr 
Ez ITS 的 测序 数据 , 利用 群体 遗传 学 和 谱系 地 理学 的 分 析 方 法 , 对 以 上 问题 进行 了 初步 的 探 
讨 和 回答 。 
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1 材料 与 方法 
1.1 材料 

本 研究 中 我 们 共 采 和 集 甘 肃 \ 四 川 、 云 南 及 贵州 27 个 居 群 的 261 个 样本 的 叶片 材料 ( 表 1; 
图 1, 图 2)， 硅胶 干燥 后 带 回 实验 室 。 

1.2 DNA 提取 、PCR 扩 增 及 测序 

依据 CTAB 法 (Doyle & Doyle, 1987) 提取 总 DNA。 对 nrDNA 的 ITS 区 段 进行 PCR 
扩 W , J 增 引物 为 ITS4(TCCTCCGCTTATTGATATGC) 和 ITS5 
(GGAAGTAAAAGTCGTAACAAGG) (White et al., 1990), X cpDNA 的 ycf6-psbM 和 
1p132-trnL 区 段 进行 PCR 扩 增 , 扩 增 引物 为 ycf6F(ATGGATATAGTAAGTCTYGCTTGGGC) 
和 psbMR ( ATGGAAGTAAATATTCTYGCATTTATTGCT ) ( Shaw et al., 2005) 及 
rpl32F(CAGTTCCAAAAAAACGTACTTC) 和 trmLR ( CTGCTTCCTAAGAGCAGCGT ) 
(Shaw et al., 2007). ITS HJ PCR 扩 增 条 件 为 95 C 变 性 2 min, 接着 以 95 C 30s, 53 °C 1 min, 
65 “C 1 min 的 条 件 循环 30 次 ， 最 后 在 65 'C 条 件 下 延伸 7 min. ycf6-psbM 的 PCR 反应 条 件 与 
ITS 基本 一 致 ， 但 其 退火 温度 为 S2 C. rp132-trnL 的 PCR 扩 增 条 件 为 80 CEHE 2 min, Z 
着 以 95 °C 1 min，50 'C Imin, 65 °C 1.5 min 的 条 件 循环 33 次 , 最 后 在 65 C 条 件 下 延伸 5 min. 
PCR 产物 用 上 海 生 工 生物 技术 有 限 公 司 的 纯化 试剂 盒 W5211 进行 纯化 , 纯化 产物 用 AB13730 W 
序 仪 进行 双向 测序 。 
1.3 数据 分 析 
测序 的 原始 数据 用 DNAStar 软件 包 的 Segman (DNAStar, Inc., Madison, USA) 进行 
HHZ, H Clustal X (Thompson et al., 1997) 进 行 比 对 。 用 PAUP*4.0b10 (Swofford, 2002) 软 
件 将 两 个 叶绿体 片段 的 序列 进行 联合 。 用 DnaSP 5.10 (Librado & Rozas, 2009) 统计 单 倍 型 
数目 和 变异 位 点 特征 ， 以 及 单 倍 型 多 样 性 (8d) 和 核 蔡 酸 多 样 性 (x); 利用 该 软件 进行 中 性 检 
验 ， 包 括 Tajima’s D 和 Fu and Li’s F* 检 验 ， 在 居 群 扩张 模型 下 推测 居 群 大 小 变化 。 运 用 
ArcGis10.2 编辑 各 单 倍 型 在 地 形 图 上 的 地 理 分 布 。 利 用 Network5.011 
(http://www.fluxus-engineering.com) 构 建 单 倍 型 网 状 进化 关系 。 

用 Permut (Pons & Petit, 1996) 软件 计算 计算 居 群 内 平均 遗传 遗传 多 样 性 Hs、 总 的 遗 
传 多 样 性 Ht 和 居 群 间 遗 传 分 化 系数 Gst 和 Nst 值 。 

应 用 Arlequin 3.5 (Excoffier & Lischer, 2010) 软件 包 中 的 分 子 变 异 分 析 (AMOVA) 分 别 
检测 该 研究 物种 在 居 群 间 和 居 群 内 的 遗传 变异 水 平 。 利用 该 软件 对 单 倍 型 分 布 的 Fst 进行 评 
价 ， 同 时 利用 该 软件 包 中 的 Mantel 统计 学 检验 ， 比 较 地 理 距 离 矩阵 与 平均 遗传 距离 窍 阵 之 
间 的 相关 性 , 并 进行 1.000 次 重复 的 显著 性 检验 ;此 外 , 失 配 分 析 (Cmismatch distribution analysis) 
也 由 该 软件 来 完成 , 即 分别 计 算 歧 点 观测 值 与 期 望 值 的 方差 (SSD)、 米 度 指数 (xr: raggedness 
index) 以 及 它们 的 显著 性 (P value)。 


表 1 样本 采集 地 信息 及 各 居 群 遗传 多 样 性 参数 统计 值 


Table 1 Information of sampling location and the values of genetic diversity index 


居 采样 地 经 度 / 纬 度 海 tk cpDNA 单 倍 型 cpDNA ER cpDNA 单 倍 ITS 单 倍 型 组 成 (序列 ITS 核 背 酸 多 态 ITS 单 倍 型 多 

群 Location Longitude Altitude (序列 数 ) 多 态 性 型 多 态 性 数 ) 性 态 性 

号 /Latitude (m) cpDNA cpDNA cpDNA ITS Haplotype cpDNA ITS 

Pop Haplotype nucleotide Haplotype (sequence number) nucleotide Haplotype 

ID (sequence diveristy (m) Diversity (Hd) diveristy (z) Diversity (Hd) 

number) 

1 甘肃 文 县 104.67°/33.05° 1 700* HI1(9) 0 0 HI) H2(2) 0.002 53 0.837 

Wenxian, Gansu H4(2) 
H15(5)* 

2 四 川 小 金 102.5°/31.02° 2 300 H1(10) 0 0 H3(10) H12(5)* 0.001 68 0.679 
Xiaojin, Sichuan H13(1)* H14(4) 

3 qJ HEY 100.99°/30.03° 2 700 H1(12) 0 0 H3(24) 0 0 
Yajiang, Sichuan 

4 四 川 乡 城 99.87°/28.75° 2 800 H1(2) H7(8) 0.000 61 0.356 H3(20) 0 0 
Xiangcheng, Sichuan 

5 UJ 99.53°/28.03° 3 180 H1(10) 0 0 H3(20) 0 0 
Derong, Sichuan 

6 云南 中 向 99.75°/28.08° 3 600” H1(4) H7(7) 0.000 88 0.509 H3(22) 0 0 
Zhongdian, Yunnan 

7 云南 中 向 100.03°/27.5° 2 630 H1(5) H7(5) 0.000 96 0.556 H3(19) H9(1) 0.000 16 0.1 
Zhongdian, Yunnan 

8 云南 玉龙 100.23°/27.29° 1 650 H1(8) H7(2) 0.000 57 0.327 H3(20) H8(2) 0.000 54 0.173 
Yulong, Yunnan 

9 云南 维 西 99.25°/27.19° 3 000” H1(10) 0 0 H3(20) 0 0 
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11 


12 


13 


14 


15 


16 


17 


18 


19 


20 


21 


Weixi, Yunnan 


云南 维 西 


Weixi, Yunnan 


云南 怒江 


Nujiang, Yunnan 


Yu 


Yo 


long, Yunnan 


云南 永 仁 


ngren, Yunnan 


BM Tie 


Ninglang, Yunnan 
中 川 盐 源 


Yanyuan, Sichuan 


Mi 


Jiu 


g/l jL 


URF 


anning, Sichuan 


long, Sichuan 


云 


和 昭通 


Zhaotong, Yunnan 


hi 


APR 


Luquan, Yunnan 


Xu 


At f 


ndian, Yunnan 


99.33°/27.12° 


98.7°/25.79° 


99.88°/26.63° 


100.4°/26.85° 


101.47°/26.42° 


100.84°/27.4° 


101.23°/27.72° 


28.4°/102.02° 


101.72°/28.53° 


103.72°/27.45° 


102.7°/26.17° 


102.84°/25.9° 


2 330 


2 700" 


2 380 


2 100 


1 400* 


2 740 


3 050 


2 110 


1 600 


2 160 


2 050 


1 100 


HI(10) 


H1(9) 


H1(2) H2(5) 


H1(8) H3(1) 


H9(1) * 


H1(8) 


H1(9) H3(1) 


H1(10) 


H1(5) H6(6) 


H1Q9) H6(1) 


H4(2) H11(8) 


H4(10) 


H7(1) H8(11) 


0.000 82 


0.000 52 


0.000 35 


0.001 04 


0.000 31 


0.000 29 


0.200 


0.167 


H3(20)) 


H3(15) H14(3) 


H1(5) H2(5) 


H2(1) H3(16) H8(2) 


H10(1) 

HIG) H2(1) 

H4(4) H5(4) 
H6(2) * 


H3(2) 


H3(16) H8(2) H9(2) 


H3(14) H8(6) + H22(2) 


H1(14) H16(5)* 


H17(3) * 


H1(4) H8(2) H18(8) * 
H19(6) * H20(4) * 


H1(9) H2(5) 
H5(5) 
H1(16) H7(4) 


H1(22) H7(2) 


H4(1) 


0.000 83 


0.001 47 


0.001 44 


0.003 03 


0.000 83 


0.001 39 


0.009 07 


0.005 85 


0.002 13 


0.000 53 


0.000 25 


0.714 


0.363 


0.858 


0.358 


0.442 


0.550 


0.737 


0.705 


0.337 


0.159 


202003.00060v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


22 云南 东 川 103.27°%/25.85° 2150 H2(1) H8(8) 0.002 06 0.378 H1(17) H7(3) 0.000 42 0.268 
Dongchuan, Yunnan H10(1) * 

23 ”云南 东 川 103.08°/25.9° 2 500 H1(1) H4(8) 0.000 19 0.222 H1(16) H11(2) 0.000 33 0.209 
Dongchuan, Yunnan 

24 ”云南 昆明 102.78°/25.32° 2 000* H1(1) H4(4) 0.002 57 0.644 H1(19) H2(1) 0.000 16 0.1 
Kunming Yunnan H5(5) * 

25 云南 澄江 102.88°/24.73° 2650 H4(10) 0 0 H1(18) H11(2) 0.000 30 0.189 
Chengjiang, Yunnan 

26 ”云南 元 江 101.8°/23.65° 1 760 H11(6) 0 0 H1 (12) 0 0 
Yuanjiang, Yunnan 

27 MÁR 108.1°/27.1° 850 H1(2) H3(3) 0.001 04 0.600 H2(2) H3(4) H8(2) 0.00391 0.8 
Shibing, Guizhou H21(2) * 

VE: * 表 示 居 群 特 有 单 倍 型 ，# 表 示 此 海拔 高 度 由 GoogleEarthPro 查询 得 到 。 


Note: *means private haplotypes; # means data estimated using GoogleEarthPro. 
p plotyp g g 


2 结果 与 分 析 

2.1cpDNA 数据 分 析 

共有 260 个 样本 的 ycf6-psbM 和 rp132-trnL 同时 测序 成 功 ， 两 个 片段 联合 并 比 对 后 矩阵 
长 度 为 1 179 bp. DnaSP 5.10 分 析 共 生成 11 个 单 倍 型 ， 每 个 居 群 的 核 苷 酸 多 样 性 (z). ë 
倍 型 多 样 性 〈Ed) 见 表 1。 物 种 水 平 单 倍 型 多 态 性 指数 Hd 为 0.656， 核 背 酸 多 态 性 (z) 为 
0.001 61. 

单 倍 型 Hl ERE ATED PG. DO) PE RBA Ae HR CA Te EE 11 个 居 
群 中 ， 昆 明 〈 居 和 群 24) 及 东 川 〈 居 群 23) 有 少量 分 布 ， 单 倍 型 H4 则 分 布 在 云南 中 部 的 5 
SR; HS 仅 分 布 在 云南 昆明 的 居 群 24 rh, Ho 仅 分 布 在 丽江 玉龙 县 的 居 群 13 中 ，H10 
仅 分 布 在 昆明 东 川 居 群 22 中 (图 1， 表 1)。Network 分 析 表 明 ， 发 生 频 率 最 高 的 单 倍 型 H1 
和 H4 位 于 网 状 结构 的 中 间 位 置 ， 单 倍 型 H3、H10、H11、H8、H5、H6 和 H5 位 于 网 状 结 
构 的 末端 位 置 (图 3)。 

居 群 的 SSD 值 以 及 糙 度 指数 > 值 都 为 不 显著 的 正 值 (CP > 0.05), 中 性 检验 表明 Tajima’s 
D (-0.59) 与 Fu and Li’s F (-0.06) 值 虽然 均 为 负 值 ， 但 不 显著 小 于 零 (P> 0.05)， 失 配 分 布 
> 曲线 为 多 峰 曲 线 〈 图 5)， 观 测 值 背离 了 期 望 值 ， 违 背 了 居 群 扩张 模型 。 这 些 结果 都 表明 异 
© 叶 泽 兰 居 群 未 发 生 显著 居 群 扩张 。 
之 异 叶 泽 兰 总 遗传 多 样 性 (Ht) 为 0.682， 居 群 内 平均 遗传 多 样 性 (Hs) 为 0.219， 居 群 
间 遗 传 分 化 系数 Gst 和 Nst 值 分 别 为 0.679 和 0.655,Nst 小 于 Gst.AMOVA 分 析 结 果 表 明 ， 
异 叶 泽 兰 较 多 的 遗传 变异 (65.55%) 发 生 在 居 群 间 , 仅 34.45% 的 变异 发 生 在 居 群 内 ,固定 指 
AX Fst= 0.655。 遗 传 距离 和 地 理 距 离 呈正 相关 ， 但 显著 度 不 高 (r= 0.18, P= 0.045)。 


2.2 ITS 分 析 
共有 261 个 样本 的 ITS 测序 成 功 ， 其 中 197 个 体 是 纯 合 体 ， 其 余 64 个 为 杂 合体 。 序 列 
w. 比 对 后 的 数据 矩阵 长 度 为 638 bp, 经 软件 DnaSP5.10 的 “phase” 功 能 拆 分 后 生成 522 条 序列 。 
= 这 522 条 序列 共生 成 21 种 基因 型 (HL-H21) (图 2， 表 1), 物种 水 平 基因 型 多 态 性 指数 Hd 
AS 为 0.686 9, FMR ASHE 0.002 35; 每 个 居 群 的 核 昔 酸 多 样 性 (z) 和 基因 型 多 样 性 (Hd) 
值 见 表 1。 
一 分 布 频率 最 高 的 为 基因 型 H1 和 H3， 其 中 HI 分 布 在 云南 中 部 及 四 川 的 晃 宁 和 九龙 地 
© 区 ， 云 南 中 部 为 其 主要 分 布地 区 ; 而 H3 从 甘肃 文 县 沿 四 川 盆 地 西部 边缘 一 直 分 布 到 分 云南 
西北 部 地 区 , 贵州 施 秉 居 群 中 也 有 分 布 , 其 中 云南 西北 部 及 四 川西 南部 为 主要 分 布 区 (图 2, 
表 1)。 基 因 型 H2 分 布 在 居 群 1、12、13、14、19、24 和 27 中 ， 除 甘肃 文 县 〈 居 群 1) 外 ， 
呈现 从 云南 西北 部 到 云南 中 部 到 云南 中 北部 ， 再 到 贵州 〈 居 群 27) 的 分 布 格局 。H4 和 H5 
呈现 出 在 甘肃 文 县 〈 居 群 1)、 云 南西 北部 〈 居 群 14) 和 云南 中 北部 〈 居 群 19) 的 分 布 式 
样 。H7 仅 分 布 在 云南 中 北部 的 三 个 居 群 中 ( 居 群 20 F121, 22), H1 仅 分 布 在 云南 中 北部 
的 两 个 居 群 中 ( 居 群 23 和 居 群 25)， 基 因 型 H8 分 布 在 四 川西 南部 ( 居 群 16)、 云 南西 北部 
( 居 群 8、13 和 15〉 以 及 贵州 地 区 ( 居 群 27)。H6 为 居 群 14 所 特有 ，H10 为 居 群 13 所 特 
有 ，H12 和 H13 仅 分 布 在 居 群 2 中 ，HI15 仅 分 布 在 居 群 1 中 ，H16 和 H17 仅 分 布 在 居 群 17 
中 ，H18、H19 和 H20 仅 分 布 在 居 群 18 中 ，H21 仅 分 布 在 居 群 27 中 。 因 此 ，ITS 基因 型 主 
要 分 布 模式 可 以 归纳 为 : A. 以 H3 为 代表 的 “甘肃 -川西 南 - 汀 西北 -贵州 ”地 区 分 布 类 型 ， 
DLAI 为 代表 的 “川西 南 - 滇 中 ” 地 区 分 布 类 型 ; B.“ 四 川 贫 地 边缘 ”特有 基因 型 分 布 类 型 ， 
如 H12, H13, H15 等 。 出 现 频率 较 高 的 单 倍 型 Hl 和 H3 位 于 网 状 进化 图 的 中 间 位 置 ， 同 
时 网 状 进化 图 出 现 多 个 环 状 结构 ， 说 明 核 基因 重组 事件 的 发 生 ( 图 4)。 
居 群 的 SSD 值 以 及 糙 度 指数 r 值 都 为 不 显著 的 正 值 (P >0.01); 中 性 检验 表明 Tajima’s 


ray 


202003.00060v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


D (-1.14) 为 负 值 ， 但 它 不 显著 小 于 零 (P > 0.1)，Fu and Li’s F 为 0.64， 失 配 分 布 曲线 为 双 峰 
曲线 (图 5S)， 观 测 值 背 离 了 期 望 值 ， 违 背 了 居 群 扩张 模型 。 由 此 推断 居 群 未 发 生 过 扩张 而 
是 处 于 动态 平衡 ， 这 与 叶绿体 DNA 的 分 析 结果 相 一 致 。 
ITS 数据 分 析 显 示 异 叶 泽 兰 总 的 遗传 多 样 性 (HL) 为 0.712, 居 群 内 平均 遗传 多 样 性 (Hs) 

为 0.325， 居 和 群 间 遗 传 分 化 系数 Gst 值 为 0.543 、Nst (HN 0.370, Nst 小 于 Gst。 AMOVA 
分 析 结 果 表 明 ， 异 叶 泽 兰 58.44% 的 遗传 变异 发 生 在 居 群 间 ， 其余 41.56% 发 生 在 居 群 内 ， 
定 指数 Fst = 0.584. EJET , 统计 分 析 显 示 异 叶 译 兰 遗传 距离 和 地 理 距 离 呈 显著 正 相 关 (Cr = 0.43, 
P=0). 
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图 1 异 叶 泽 兰 叶绿体 单 倍 型 地 理 分 布 。 
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Fig. 1 Geographical distribution of combined cpDNA haplotypes of Eupatorium heterophyllum 
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图 2 异 叶 泽 兰 ITS 基因 型 地 理 分 布 。 
Fig.2 Geographical distribution of ITS genotypes of Eupatorium heterophyllum 
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圈 代 表 一 种 单 倍 型 ， 圆 圈 的 大 小 与 单 倍 型 的 频率 成 正比 ， 每 一 分 支 代 表 一 次 突变 ， 黑 色 的 
小 圆圈 代表 缺失 单 倍 型 (mv)。 下 同 。 


Note: Each circle represents one haplotype, the size of circles corresponds to the frequency of each haplotype, each 


$: 图 中 每 个 


hay 


branch represents one mutation, and the black open circle represents hypothetical missing haplotype (mv). The 
same below. 
3 叶绿体 DNA 联合 分 析 的 单 倍 型 网 状 进化 树 
Fig.3 Network analysis of the combined cpDNA sequences 
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Fig.4 Network analysis of the ITS sequences 
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5 叶绿体 DNA 序列 数据 的 岐 点 分 布 〈 左 ) 和 ITS 序列 数据 的 岐 点 分 布 〈 右 ) 
Fig.3 Mismatch distribution for cpDNA data (left) and mismatch distribution for ITS data (right) 
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3 讨论 
3.1 异 叶 泽 兰 的 遗传 多 样 性 和 遗传 结构 
本 研究 中 ， 我 们 对 27 个 居 群 的 261 个 个 体 的 三 个 cpDNA 片段 及 ITS 进行 了 测序 和 分 
析 ， 三 个 cpDNA 片段 联合 分 析 后 共生 成 11 个 单 倍 型 ，ITS 序列 经 PHASE 拆 分 后 共生 成 21 
个 单 倍 型 ， 表 现 出 较 高 的 单 倍 型 多 样 性 。 在 物种 水 平 上 ，cpDNA BJE BR AVE (z) 为 
0.001 61, ITS 序列 核 苷 酸 多 态 性 为 0.002 3$。 在 菊 科 植物 中 ， 内 蒙古 革 苞 菊 Tugarinovia 
mongolica 是 菊 科 一 个 单 种 属 植物 ,包含 一 变种 Tugarinovia mongolica var. ovatifolia,， 原 变种 
分 布 在 内 蒙古 北部 地 区 ， 其 卵 叶 变 种 分 布 在 内 蒙古 南部 。 基 于 两 个 叶绿体 DNA 片段 
(psbA-trnH 和 psbK-psb1) 的 数据 分 析 表 明 ， 内 蒙古 革 苞 菊 具 有 非常 高 的 核 昔 酸 多 样 性 (x 
= 0.009 2) 和 单 倍 型 多 样 性 (Ed = 0.908 6)， 两 个 变种 之 间 遗 传 分 化 明显 (Zhao et al., 2019). 
白菊 木 (Leucomeris decora) 为 菊 科 落叶 小 乔木 ， 其 trnO-rps16, rp116 和 rp132-trnL 联合 数 
HR BJ ARS ASPEN z = 0.001 02， 而 变异 率 相对 较 高 的 核 基因 片段 GAPDH ËJ z = 0.002 37 
(Zhao & Gong, 2012). RSE (Ligularia hodgsonii) 为 中 国 - 日 本 间断 分 布 的 菊 科 窜 吾 属 
植物 ,在 中 国 主要 是 围绕 四 川 盆 地 分 布 , 其 23 个 居 群 的 三 个 叶绿体 片段 联合 分 析 
(trnQ-5'rps16, trnL-rpl32 和 psbA-trnH) WKF AR AS VE zx = 0.002 99， 单 倍 型 多 态 性 为 Hd = 
0.847(Wang et al., 2013)。 因 此 ， 与 这 些 菊 科 植物 相 比 ， 异 时 泽 兰 物种 水 平 的 核 昔 酸 多 态 性 和 
单 倍 型 多 态 性 都 相对 较 低 。 但 是 ， 在 居 群 水 平 ，ITS 数据 显示 四 川 盆 地 边缘 的 几 个 居 群 存在 
相对 较 高 的 单 倍 型 多 样 性 ( 表 1， 图 2)。 
cpDNA 单 倍 型 数据 显示 居 群 水 平 总 的 遗传 多 样 性 (Ht = 0.764) 大 于 居 群 内 平均 遗传 多 
样 性 (Hs = 0.250)， TIRE, ITS 数据 显示 总 的 遗传 多 样 性 (Ht = 0.746) 也 大 于 居 群 内 平均 遗 
传 多 样 性 (Hs = 0.382). AMOVA 分 析 表 明 , 异 叶 泽 兰 居 群 间 遗 传 变异 高 于 居 群 内 (cpDNA: 
Fst = 0.655; nITS: Fst = 0.584)。 而 Permut 分 析 显 示 居 群 水 平 的 遗传 分 化 系数 相对 较 高 
(cpDNA: Gst = 0.679, Nst= 0.655; nITS: Gst = 0.543, Nst = 0.370)。Petit et al. (2005) 
曾 统计 了 124 种 被 子 植物 的 属于 母性 遗传 的 分 子 标 记 遗 传 分 化 系数 Gst 平均 值 和 77 种 被 子 
植物 的 属于 双亲 遗传 的 分 子 标记 遗传 分 化 系数 Gst 平均 值 。 结 果 显 示 ， 属 于 母性 遗传 的 遗传 
分 化 系数 Gst 平均 值 为 0.637， 而 属于 双亲 遗传 的 遗传 分 化 系数 Gst 值 平均 为 0.184， 异 叶 泽 
兰 分 化 系数 高 于 这 两 个 值 . 以 上 数据 都 表明 , 异 叶 泽 兰 居 群 间 存 在 较 明 显 的 遗传 分 化 。 然 而 ， 
由 于 Ast 值 小 于 Gst 值 ， 异 叶 泽 兰 的 分 布 不 具有 明显 的 谱系 地 理 结构 (Pons & Petit, 1996). 
但 是 ， 从 单 倍 型 分 布 式样 看 ， 却 存在 着 较 明 显 的 “甘肃 -川西 南 - 汗 西 北 -贵州 ”和 “ 污 中 ” 
两 个 相对 分 离 的 分 布 区 ,两 个 分 布 区 都 存在 各 自 特 有 单 倍 型 和 主要 单 倍 型 。 然 而 单 倍 型 网 状 
进化 树 显 示 , 两 个 分 布 区 所 拥有 的 单 倍 型 并 没有 各 自 形成 两 大 谱系 分 支 , 同时 两 个 分 布 区 间 
也 存在 分 布 频率 相对 较 低 的 共享 单 倍 型 。 作 为 一 种 多 年 生 菊 科 植物 , 异 叶 泽 兰 具有 较 强 的 适 
应 能 力 ,而 种 子 具有 较 强 的 远 距 离 传播 能 力 ， 两 分 布 区 间 应 该 存在 较 强 的 基因 交流 ,这 可 能 
是 导致 异 叶 泽 兰 遗 传 多 样 性 分 布 式 样 的 主要 原因 。 
3.2 异 叶 泽 兰 的 居 群 历史 动态 
谱系 地 理学 研究 的 主要 目的 之 一 是 推测 第 四 纪 冰 期 时 某 物 种 的 避难 所 及 冰 后 期 分 布 范 
围 的 扩展 (Avise, 2000; Liu et al., 2012)。 青 藏 高 原 东 部 及 邻近 的 中 国 西南 地 区 ， 分布 着 大 量 
的 古老 种 和 新 起 源 种 ， 其 中 许多 植物 的 分 布 范 围 在 冰 后 期 发 生 过 扩张 (Liu et al., 2012; Qiu 
et al., 2011), 但 也 有 未 生 扩 张 的 物种 , 如 偏 花 报 春 (Primula secundiflora) Wang et al., 2008). 
而 异 叶 泽 兰 作为 分 布 在 青藏 高 原 东部 和 横断 山脉 的 物种 , 中 性 检验 和 失 配 分 析 表 明 其 未 发 生 
过 居 群 的 扩张 。 
冰期 时 的 避难 所 往往 会 保留 较 高 的 遗传 多 样 性 和 单 倍 型 多 样 性 ， 而 单 倍 型 分 布 能 够 反 
映 一 个 物种 的 地 理 分 布 格局 特征 (Hewitt, 1996, 2000;Petit et al., 2003)。 根 据 漳 祖 理论 ， 分 布 


频率 较 高 的 , 位 于 网 状 进 化 树 中 间 位 置 的 单 倍 型 可 能 为 较 古 老 的 单 倍 型 ， 而 居 群 特有 单 倍 型 
可 能 是 通过 近期 辐射 分 化 衍生 而 来 的 较 年 轻 的 单 倍 型 (Emerson et al., 2001)。 依 据 Network 
分 析 结 果 ， 异 叶 泽 兰 叶绿体 DNA 两 个 单 倍 型 Hl 和 H4 以 及 ITS 两 个 单 倍 型 Hl 和 H3 分 布 
的 地 理 范 围 较 广 ， 且 位 于 网 状 结构 的 内 部 位 置 ， 可 能 为 古老 单 倍 型 (Templeton et al., 1992; 
Crandall & Templeton, 1993)。 其 中 ， 叶 绿 体 DNA 单 倍 型 Hl 和 ITS 的 单 倍 型 H3 的 广泛 分 
布 于 甘肃 -川西 南 - 滇 西北 一 带 及 贵州 地 区 ， 而 叶绿体 DNA 单 倍 型 H4 和 ITS 的 单 倍 型 H1 
则 分 布 于 云南 中 部 地 区 。 由 此 推测 ， 横 断 山 区 南部 《川西 南 - 滇 西北， 和 云南 中 部 可 能 是 异 
叶 泽 兰 在 第 四 纪 冰 期 尤其 是 末次 冰 盛 期 (LGM) 时 的 两 个 避难 所 所 在 地 。 
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